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1 Introduction 



The asymmetric simple exclusion process (ASEP) |TJI2] has often been called 
the Ising model of nonequilibrium statistical physics. In this stochastic lattice 
gas model particles move randomly with a bias onto neighboring lattice sites, 
provided the target site is empty. Even its most-studied one-dimensional 
version which describes driven single-file diffusion exhibits rich phenomena, 
in particular boundary- induced phase transitions El 03 EH Ej, and has a 
wide range of applications in different branches of physics. Experiments 
probing single-file diffusion have been performed with molecules in zeolites 
[S], colloidal particles in confined geometry |Hj or optical lattices [TU], and 
very recently with granular systems jTT] . Driven single-file transport has been 
observed in biological systems and the ASEP serves as a starting point for 
modelling the motion of ribosomes along the m-RNA during protein synthesis 
[T2*l ITS] and of molecular motors along microtubuli or actin filaments ^3] . 
The ASEP is also at the heart of the cellular automaton approach to vehicular 
traffic on road networks ^Hl EH] • In this setting some predictions from the 
theory of boundary-induced phase transitions for the development of traffic 
jams have been verified empirically despite the complexity of real vehicular 
traffic ^7j. Phenomena analogous to traffic jams have also been observed in 
the biological transport systems referred to above. 

In the hydrodynamic approach to traffic flow 16 using partial differential 
equations traffic jams correspond to shock solutions. A shock is a density 
discontinuity on moving with some deterministic speed, determined by mass 
conservation. It is therefore no surprise that on macroscopic Euler scale the 
time evolution of the particle density of the ASEP is described by the inviscid 
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Burgers equation [TB1HHI which develops shocks for generic initial data. With 
a view on applications of the ASEP to systems for which a hydrodynamic 
description is too coarse-grained it would thus be of interest to understand 
what fluctuating microscopic structure (on lattice scale) is underlying the 
phenomenon of shocks. 

In fact, a great deal is known about shocks in the ASEP due to the exact 
solubility of the model. In the stationary regime the shock structure has 
been studied as seen from a so-called second-class particle which serves as 
microscopic marker of the shock position. The particle density decays away 
from the shock exponentially (on lattice scale) to the respective constant 
bulk values pip of the two branches of the shock [201 EH 122] ■ The shock 
position itself has been proved to perform Brownian motion on coarse grained 
diffusive scale j2Bj- For a particular strength of the driving field the associated 
decay constant of the particle density vanishes, corresponding to a "minimal" 
intrinsic shock width. For this special value of the driving field also the 
motion of the shock simplifies greatly. It performs a biased random walk on 
microscopic lattice scale with explicitly known hopping rates I25j. 

It is natural to ask whether this special feature of the ASEP survives in 
more complicated models of driven diffusive systems. In particular, one would 
like to investigate exclusion processes with nonconserved internal degrees of 
freedom, where particles may have different velocities, charges, masses or 
other distinguishing properties that they can gain or lose e.g. in a collision 
or chemical reaction. Here we address this question in the simplest case 
of two possible internal states that each particle can posses. Such models 
have been investigated recently for various biological and vehicular transport 
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phenomena f3J I2S1 I2Z1 • Studying the microstructure of a shock illuminates 
the role of finite-size effects in first-order boundary-induced phase transi- 
tions that are associated with the motion of traffic jams [2H1 EDI ED! in finite 
systems. 

The paper is organized as follows: In the following section we determine 
the families of three-states models with pair interaction and a single con- 
servation law. We also define shock measures for these systems. In Sec. 3 
we study exclusion processes with binary internal degree of freedom that al- 
lows for special travelling shock solutions on the finite lattice. In Sec. 4 we 
summarize our results and draw some conclusions. 

2 Three-states processes with one conserva- 
tion law 

2.1 Stochastic dynamics 

On an abstract level the exclusion process with a binary internal degree of 
freedom is a three-states process where the state of the system at any given 
time is described by a set of "occupation numbers" n = m, . . .,ni where 
rik = 0, 1, 2 is the local occupation number at site k and L is the number of 
sites. In the next section we assign state to an empty lattice site, state 1 to 
a particle of type A and state 2 to a particle of type B. The labels A and B 
represent two possible internal states of a particle. However, in this section 
we first keep the treatment general and consider the occupation numbers as 
abstract objects, labelling one out of three possible states of a lattice site. 

The bulk stochastic dynamics are defined by nearest neighbor transitions 
between the occupation variables which occur independently and randomly 
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in continuous time after an exponentially distributed waiting time. The 
mean r(n' k , n' k+1 ; rik+i) of this waiting time depends on the transition 
(rik,nk+i) — > (n' k ,n k+1 ), its inverse is the transition rate. 

At the boundary sites 1, L we assume the system to be connected to 
some external reservoir with which the system can exchange particles. For 
injection and extraction of particles at the left boundary we introduce the 
rates : 

0^1 «i, 71, 
0^2 a 2 , 72, 

1-2 a 3 , 73, (1) 

and for the right boundary 

0^1 6 U /3 U 
0^ 2 5 2 , fo, 

1 # 2 5 3 , /3s. (2) 

Here and below the left rate refers to the process going from left to right, 
while the right rate is for the reversed process. 

The time evolution defined above can be written in terms of a continuous- 
time master equation for the probability vector 

\P(t))=J2 p (mr--,n L ;t)\n), (3) 

n 

where P(rii, • ■ ■ , n^; t) is the distribution for the probability of finding parti- 
cles at sites 1 to L and | n) is the basis vector in the space of configurations 
in the naturally defined tensor basis |2j. The probability vector is normalized 



such that (s\P) = 1 with the summation vector (s\ — I] n (n| and scalar 
product (n\n') = Sn,n'- The time evolution is generated by the stochastic 
"quantum Hamiltonian" H whose offdiagonal matrix elements H^n' are the 
negative transition rates between configurations. As required by conserva- 
tion of probability, the diagonal elements are the negative sum of transition 
rates in the respective column. 

Therefore the master equation is now described by the imaginary time 
Schrodinger equation 

j t \P(t)) = -H\P(t)) (4) 

with the formal solution 

|P(f)>=e-*'|P(0)>. (5) 

Since only nearest-neighbour interactions are included, the quantum Hamil- 
tonian H defined above has the structure 

L-l 

H = b 1 + Y. Kk+i + h- (6) 
fc=i 

where b\ and are the boundary matrices 

+ a 2 ) 7i 72 
6i = - [ an -(71 + a 3 ) 73 j . (7) 

«2 «3 -(72 +73) 



b L = - S, -(Pi + 6 3 ) . (8) 

V $2 P3 -(P2+Pl)) L 

The local bulk transition matrix hk,k+i acts non-trivially only on sites k 
and k + 1. To define its matrix elements we introduce an integer label 



3n k + n k+1 + 1 



(9) 
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in the range 1 < i < 9 for the occupation variables on two neighboring 
sites. The offdiagonal matrix elements (hk,k+i) are then the transition 
rates —w%j. Here i = 3n' k + n' k+1 + 1 labels the target configuration and j is 
the respective label of the initial configuration (nfc,n fe+1 ). 

2.2 Symmetries and conservation laws 

Within this setting one could describe 72 different bulk transitions, corre- 
sponding to the 72 mathematically possible changes of configurations on a 
pair of sites. However, we shall reduce this large number by imposing various 
physically motivated constraints. First, we require a local conservation law. 
Generally, the physical interpretation of the conservation law depends on the 
physical interpretation of the occupation numbers nu and will become clear 
below. Mathematically this means that in a periodic system some function 
J2k C(nk) of the local occupation numbers should remain invariant under the 
stochastic dynamics, i.e., 

CK) + C« +1 ) = C(n k ) + C(n k+1 ) (10) 

for any local transition between configurations This constraint forces 
a large number of transition rates tUy to vanish. Physically C(n) is some 
observable property (such as mass or charge) of the state n. 

The conservation condition (jl(J|) does not uniquely define the function 
C(n). In order to analyse these constraints we set C(0) = and C(l) = 1. 
This involves no loss of generality since adding a constant to C(n) or mul- 
tiplying C(n) by an arbitrary factor leaves (jl(J|) invariant. From inspec- 
tion of (jl(Jj) one can then see that there are three distinct families of solu- 
tions: (i) degenerate case, represented by C(2) = C(l) = 1 (or equivalently 
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C(2) = C(0) = 0), (ii) linear nondegenerate case, represented by (7(2) =2 
(or equivalently C(2) = —1, C(2) = 1/2), (iii) two independent conservation 
laws, represented by any other value of C(2). The non-degenerate linear con- 
servation law is treated elsewhere (231, the case of two conservation laws was 
studied in detail in [3TJ E2] . Here we investigate the degenerate conservation 
law. The degenerate function C(n) has a natural interpretation as counting 
the number of particles at a given site irrespective of its internal state. This 
is the motivation behind the assignment of the state labels A, B used below. 
The presence of a conservation law implies a lattice continuity equation 

^C fc =j fc -i-jfc (11) 

for the expectation Ck = (C(n^)). This quantity plays the role of a local 
order parameter. The quantity jk is the current associated with the con- 
servation law. It is given by the expectation of some combination of local 
occupation numbers, depending on the model under investigation, see below. 
Since we do not study here periodic systems we do not require the boundary 
sites where the system is connected to the reservoir to respect the conserva- 
tion law. The quantities jo, ji entering the continuity equation for k = 1 
and k = L respectively are source terms resulting from the reservoirs. They 
are functions of the reservoir densities. 

Second, in addition to the conservation law we require PT-invariance, i.e., 
the bulk dynamics should be symmetric under combined time reversal and 
space reflection. This physical input generalizes the equilibrium condition of 
detailed balance to allow for external driving forces which lead to a bias in 
the hopping rates. In such the system is forced into a nonequilibrium 

steady state with a stationary current flowing in the system. Well-known 
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examples for models of this kind are exclusion processes satisfying pairwise 
balance PT- invariance is implemented by demanding detailed balance 
with respect to the reflected target state of a local transition. As a result, 
there are pairwise relations between some of the 72 transition rates, see below. 

2.3 Product measures 

Even though the number of independent model parameters is greatly reduced 
by particle conservation and PT-symmetry the form of the stationary distri- 
bution is not determined by these constraints. In order to able to carry out 
explicit computations we restrict ourselves to systems such that the station- 
ary distribution of the stochastic dynamics factorizes, i.e., one has a product 
measure. In the quantum Hamiltonian formalism introduced above a product 
measure is given by a tensor product 



Here the three-component single-site probability vectors \Pk) has as compo- 
nents the probabilities P(rik) of finding state n at site k. In the stationary 
distribution these probabilities are position- independent, \P^) = \P), and the 
stationary probability vector thus has the homogeneous product form 




(12) 



(13) 



We represent the single-site basis vectors for this model as 




(14) 
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and parametrize the stationary one-site marginal 

(15) 




by a fugacity z and the ratio c of A and B concentrations. The normalization 
factor 

v = 1 + z + cz (16) 
is the local partition function. Thus one has for this grandcanonical ensemble 

P A = ~, P B = c- (17) 

V V 

and for the total conserved particle density 

p : = / + p z = z l- In v = (1 + c )-. (18) 

In formal analogy to systems in thermal equilibrium we shall refer to the 
logarithm of the fugacity as chemical potential. 

By definition of stationarity the stationary probability vector satisfies the 
eigenvalue equation 

H\P*) = 0. (19) 

Requiring the existence of a stationary product measure imposes constraints 
both on the bulk rates and on the boundary rates which fix the bulk fugacity 
z. Once these conditions are determined the model is fully defined and 
its stationary distribution is given. Notice that by definition a stationary 
product measure has no correlations between the occupation numbers at 
different sites. 

After defining the model in this way we shall relax some of the constraints 
on the boundary conditions and study the time evolution of (nonstationary) 
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shock measures of the form 

\k) — \Pt) m ® \P 2 )® L ~ k . (20) 

These shock measures have single-site probabilities given by a fugacity z\ in 
the left chain segment up to site k and fugacity z 2 in the remaining chain 
segment from site L — k up to site L. Such a shock measure defines fully the 
internal structure of the shock. Since there are no correlations in a shock mea- 
sure one may regard the lattice unit as the intrinsic shock width. A typical 
configuration has a sharp decrease of the mean interparticle distance across 
the lattice point k. The boundary fugacities of the system are chosen such 
that each chain segment is stationary at its boundary. The measure itself, 
however, is not stationary for z\ ^ z 2 . The associated gradient of the chemi- 
cal potential together with external driving forces entering the bulk hopping 
rates drive the system into an nonequilibrium steady state, to be determined 
below as the final stage of the time evolution of the shock measure. 

3 Exclusion process with binary internal de- 
gree of freedom 

We now implement the constraints discussed above. The degenerate con- 
servation law ()10|) forces 48 transition rates to vanish. The following 24 
transitions remain: 

OA -> AO w 42 , AO — ► OA w 24 , 
OB -»■ BO w 73 , BO -> OB w 37 , 
AB -> BA w 8 e, BA -> AB w 68 , 
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BO -» 


■ AO 


w 47 , 


OA- 


> OB 




OB -» 


■ AO 




OA- 


> 50 




AO -» 


BO 




OB - 


* OA 




BO -> 


■ OA 


^2/ 1 


AO - 


> OB 




£A - 


♦ AA 


^58, 


AA - 


> A£ 


^65, 


AS - 


> AA 


W 56 , 


AA - 


> 5A 


^85, 


fifi - 


* AA 


^59, 


AA - 


> Bfi 




fifi - 


> 5A 


W 89 , 


AB ~ 


* fifi 


^96, 


BB - 


> A5 


W69, 


BA - 




w 98 . 



(21) 

Parity-time invariance leads to pairwise relations between some of these 
rates. Time reversal symmetry means to have detailed balance p*{n)w(n — > 
n') = p*(n/)w(n' n). In order to combine this relation with the parity 
(space reflection) operation we change the position of neighbouring sites with 
each other in the initial configuration and final configuration on the left-hand 
side of the detailed-balance relation. Using (|TK|) this yields the following 
symbolic relations for the rates 

w(A -»• B) = cw(B -> A) (22) 

for each particle on a pair of neighboring sites. With this relation we can 
reduce the number of independent rates in the process (|21|) to only 15 non- 
stationary rates, viz. 6 hopping rates and 9 "reaction rates" for changes of 
the internal states of the particles. For clarity we represent all of the hopping 
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rates by fa's and reaction process by r's and write the rates as 

W47 = ri, w 32 = cri, 
w 43 = r 2 , w 72 = cr 2 , 
w 23 = r 3 , w 7A = cr 3 , 
w 27 = r 4 , w u = cr 4 , 
w 58 = r 5 , w 65 = cr 5 , 
w 5 6 = r 6 , w 85 = cr 6 , 
w 59 = r 7 , w 95 = c 2 r 7 , 
w 89 = r 8 , w 96 = cr 8 , 
w 69 = r 9 , w 98 = cr 9 , 
w A2 = h u w 73 = h 2 , 
w 24 = h 3 , W S 6 = h A , 
w 68 = h 5 , w 37 = h 6 . 



(23) 



In the quantum Hamiltonian formalism, the bulk transition matrix is then 
j;iven by 



hk,k+i — 



( ■ 
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r 3 


h 3 








r 4 















cr^ 








h 6 













r 2 










7-1 






















r 6 





r 5 


r 7 














cr 5 







h 5 


r 9 





cr 2 


h 2 


cr 3 




























cr 6 


/i 4 







^8 


V 











c 2 r 7 


cr 8 





cr 9 





(24) 



k,k+l 
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3.1 Product measure 

With (|15J1 the homogeneous product measure has the form 

1 



cz 



It is convenient to define 

h i>i+1 = ht- [E(n? - h? +1 ) + E'(nf - nf +1 )] (26) 

where E, E' are arbitrary constants and ua and ub are number operators with 
eigenvalue 1 if a particle of the respective species is present and otherwise. 
Furthermore we define modified boundary matrices 

b 1 = b 1 + Enf + Enf, b L = b L - Ehf - E'nf. (27) 

This allows us to rewrite the quantum Hamiltonian as 

L-l 

+ + (28) 

»=i 

The eigenvalue equation (|19|) may be rewritten 

= h hl+1 \P*) = (b! + g)\P*) = (b L -g)\P*). (29) 

with a further arbitrary constant g. 

This trick allows us to determined the conditions on the rates that ensure 
that actually is stationary. For the bulk rates (|2*9"|) yields 

E = h 3 - hi + c(r 3 + r 4 - n - r 2 ), 

E' = h 6 - h 2 + n + r 4 - r 2 - r 3 . (30) 
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Furthermore, some algebra shows that the bulk rates must satisfy the follow- 
ing condition for stationarity 

h6-h 2 +h 1 -h 3 +h A -h 5 + (l+c)(r 1 -r 3 ) + (l-c)(r 4: -r2)+c(r 8 -r 9 )+r 6 -r 5 = 0. 

(31) 

In order to satisfy the eigenvalue equation at the boundaries the terms in- 
volving E, E' must cancel. For the left boundary this yields the two relations 

[hi - h 3 + c(2(r 2 - r 4 ) + h 2 - h 6 )]z = ^zu + ^ 2 czv - («i + a 2 )u 

= -fiizv - f3 2 czv + (5 1 + 5 2 )u, 

(32) 

and similarly at the right boundary 

[(-r 5 + r 6 + c(r 8 - r 9 ) + h 4 - h 5 )z - r l + r 2 + r 3 - r 4 + h 2 - h 6 ]cz 

= (l2 + lz)czv - a 2 v - a 3 z(l + z)v 
= -(/3 2 + (3 3 )czu + 5 2 u + 5 3 z(l +z)u. (33) 

These relations define a model for which the product measure with constant 
fugacity z is stationary. The fugacity is determined by its boundary value 
encoded in the boundary rates. 

3.2 Fugacity gradient 

Now we generalize the model to allow for different fugacities z±, z 2 at the two 
boundaries. The product measure is then no longer stationary and there is no 
general principle that would constrain the form of the stationary distribution. 
However, in principle its properties can be calculated from the studying the 
time evolution of the system starting from some initial distribution. 
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In general, solving for the dynamics of a many-particle system is a much 
harder task than determining its stationary distribution. However, guided 
by previous experience |2S] we make as ansatz an initial distribution which is 
a shock measure connecting the two boundary fugacities. The representation 
of the shock measure here is 

(g>L— k 



k) 









( 1 


k) 


(•') 


• 


Z 2 








V cz 2 



(34) 



On a coarse-grained scale the density profile corresponding to this measure 
has a jump discontinuity, see Fig. ^ We search for conditions on the rates 
such that 



d_ 

If 



k) = di\ k - 1) + d 2 \ k + 1) - (di + d 2 )\ k). 



(35) 



This implies that the family of shock measures labelled by the shock position 
k is closed under the time evolution of the many-particle system. Physically 
this behaviour corresponds to a random walk of the shock with hopping rates 
di, (d 2 ) to the left (right). 



(2) 

Pb 

(2) 

Pa 



1 k L 

Figure 1: Coarse grained density profiles of a shock measure with shock 
between sites k, k + 1. 



In order to have the random walk equation (J3*5j) for the shock, one replaces 
the left hand side by the (negative) quantum Hamiltonian in the form 
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Then in each branch of the shock one has k) = 0, except for i = k. 

Stationarity at the boundaries implies 

h\P*) = (-Ehf - E'hf +gi)\P*), (36) 

b L \P*) = {Eni + E>h B L - g 2 )\P*), (37) 
where gx and g 2 are obtained using (jBTJj) as 

91 = E^ + E'^ = (l + c)(p-q)^ 

= a x + a 2 - (71 + 072)21, (38) 

^ 2 = E^ + Ef^ = {l + c){p-q)^ 

V 2 V 2 V2 

= -(Si + 5 2 ) + (Pi + cP 2 )z 2 . (39) 

The random walk condition for the shock thus leads to 9 equations 

{-tik + dt + dz -gi + g 2 )\k) - 41 k - 1) - d 2 \ k + 1 ) = (40) 

for the bulk rates. Three of these conditions are fulfilled without any con- 
straint on the rates, leaving the following 6 equations : 

S-d 1 —-d 2 — = 0, (41) 
v 2 v x 

(cr A + hz)(zi - z 2 ) + Sz 2 - d x z 2 — - d 2 z x — = 0, (42) 

v 2 Vx 

(r 4 + h 6 )(zi - z 2 ) + Sz 2 - diz 2 — - d 2 z 1 — = 0, (43) 

v 2 v x 

(cr 2 + hi)(z 2 - z x ) + Szi - d x z 2 ^- - d 2 z 1 ^ L = 0, (44) 
-4— = 0, " ' (45) 

ZlV 2 Z 2 V X 

(r 2 + h 2 )(z 2 - z x ) + Szx - d x z 2 — - d 2 z x — = 0, (46) 

v 2 v x 
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where 

S = d 1 + d 2 + g 2 -g 1 . (47) 

Solving the above equations leads to three independent relations between 
bulk rates and densities 

h 3 + cr 4 = h 6 + r 4 = p, (48) 
hi + cr 2 = h 2 + r 2 = q, (49) 

P - = "A (50) 
q zi 

and two relations 

di = q~, (51) 

d 2 =p-, (52) 
v 2 

that express the shock hopping rates in terms of the hopping rates of the 
model and the fugacities of the shock. On this parameter manifold the sta- 
tionarity condition (pTTj) reduces to 

h± - hs + (1 + c)(ri - r 3 ) + c(r 8 - r 9 ) + r 6 - r 5 = 0. (53) 

The shock performs a random walk for a specific ratio of the boundary 
fugacities, or, equivalently, at some specific strength of the driving force 
encoded in the particle hopping rates. Thus shock mean velosity v s in terms 
of vacancy density and hopping rates is 

«. = SL^l, (54) 
v x v 2 

and its diffusion coefficient as long as the shock is far from the boundaries is 

D. = 4±5!i. (55) 
18 



From the shock hopping rates and its biased random walk dynamics we 
can read off the stationary distribution of the system for different boundary 
densities. This is a linear combination of shock measures 



For d\ > g?2 (bias to the right) the stationary shock position is in the vicinity 
of the right boundary, leaving the system in a phase of low density. Con- 
versely, for d\ < g?2, the system is in a high-density phase. At d\ = <i 2 the 
system undergoes a first-order nonequilibrium transition j^j. Here the shock 
has no bias and can be found with equal probability anywhere on the lattice. 
The stationary density profile is linear, but a typical particle configuration 
has two different regions of constant (but fluctuating) density. The density 
jumps quickly from one density to another in some small region of the lattice. 

3.3 Steady state current 

In order to make contact with the ASEP we calculate the stationary current 
for this model. In order to identify the current we first calculate the equation 
of motion for the expected local particle densities, 



-<n£) = -(/ii+cr 1 + cr 2 )(n°_ 1 n^) + / i3 (^_ 1 n°) + (/ i 4-r 6 )(^„ 1 nf) 
-(h 5 + cr^in^ni) + ^(n^nf ) + ^(nf-Xt) 
-c(r 5 + cr^ntA) + (r 7 + r 8 )«_ x nf ) + M™k<-i> 
-(h 3 + cr 3 + cr 4 )(nfn° k+1 ) - [h A + cr 8 )(nfn^ +1 ) 
+{h 5 + r 5 )(nf^ +1 ) + n{n^nl +l ) - r 2 {n k n* +1 ) 
-c(r 6 + cr 7 )(n£n£ +1 ) + (r 7 + r 9 )«nf +1 ), 




(56) 



d 



(57) 
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j t (n*) = Hh2 + r 2 + r 3 )(nl_ 1 n*) + h e (ntA)-(h 4 + r 6 )(nt 1 nf) 
+(h 5 + cr 9 )(nf_ 1 n^) + cr^n^nf) + cr^n^nf) 
+c(r 5 + cr 7 ){ntA) - (r 7 + r^n^nf) + M^f+i) 

+ n + ^Xnfnfc+i) + {hi + cr 8 )(n^nf +1 ) 
-(/is + r 5 )(n^n£ +1 ) + cr 2 (n° k n£ +1 ) + cr 3 (n£n° k+1 ) 
+c(r 6 + cr 7 ){n£n£ +1 ) - (r 7 + r 9 )(ra£n£ +1 ). 



This can be written in terms of A and S particle current 



where the source term 



'^(ntA) - (cr 5 + c 2 r 7 ){ntA) + r^ntX) 
+ (r 7 + r 8 )(ntX) ~ cr^n^n*) + (n + y)<« +1 > 
+y(«i> - - (cr 3 + ^)<nj*n° +1 > 

+r 5 (nf n£ +1 ) - (cr 6 + c 2 r 7 )(n^ +1 ) + (r 7 + r 9 )(nf nf +1 ) 
-cr 8 (ra£n£ +1 ). 



(58) 



l(ni)=j£-i-tf + Sk, (59) 
^(nZ)=j?-i-ti-Sk, (60) 



(61) 



expresses the fact that the individual particle densities are not conserved. 
The particle currents are given by the expectations 

3k = Hhi + ^)(nlni +1 ) + (h 3 + ^)(nin° k+1 } + h 4 (nin^ +1 } 
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jj? = -(^2 + ^)(« +1 ) + (/i 6 + ^)(nfn° +1 )-/i 4 (n^nf +1 ) 
+hs(n%n£ hl ) - C ^(nlni +l ) + ^i( n ^n° +1 ). 



(62) 



(63) 



By adding the two individual currents we find the total particle current 
to be given by 

3 k — 3 k + 3 k 

= -/iiK<-i> + h 3 (n£n° k+1 ) - h 2 (n° k n* +1 ) + h 6 (n*n° k+1 ) 
-r 2 (n° k n^ +1 ) - cr 2 (n° k nf +1 ) + cr 4 (nfn° k+1 ) + n(nfn° k+1 ). 

(64) 

In the steady state we obtain 

h 3 -h 1 + c(h 6 - h 2 ) + 2c(r 4 - r 2 ) 

3 = 7— P(l-P), (65) 

1 + c 

where p is the average density (JTHJ). This can be written in terms of E and 
E' 

, E + cE' 

f = -±^-p(l-p). (66) 

This is the well-known parabolic current-density relation of the ASEP PjJ |2] 
where the density-independent prefactor plays the role of the hopping bias. 
In fact, on the special manifold which gives rise to the random walk of the 
shock we find, using (pf%j) - (}4U|) . the simpler expression 

j* = (p-q)p(l-p). (67) 
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4 Conclusions 



We have found that three-states lattice gases with a single local conservation 
law can be classified into two families, one where the conserved quantity 
is a linear function of the occupation variable, another where the function 
is degenerate, i.e., takes the same value for two different states. Nonlinear 
nondegenerate functions lead to two independently conserved quantities. 

The degenerate linear conservation describes a class of asymmetric exclu- 
sion processes with a binary internal degree of freedom. We have identified 
constraints on the transition rates such that the stationary distribution is a 
product measure, parametrized by the nonequilibrium analog of the fugac- 
ity For open systems with different boundary fugacities we have found a 
complete list of models where the shock performs a biased random walk on 
the lattice. For these systems we have detailed knowledge about the micro- 
scopic structure of the shock. As in other models studied previously (see [33] 
and references therein) these shocks are intrinsically maximally sharp and 
behave like collective single-particle excitations already on the lattice scale - 
not only after coarse-graining where all the microscopic features of the shock 
are lost. Apparently this enormous reduction in the number of dynamical 
degrees of freedom in a subspace of the stochastic dynamics appears more 
frequently than previously suggested [35]. 

An immediate consequence of the random walk dynamics of the shock is 
the existence of a first order boundary- induced phase transition which occurs 
if the boundary fugacities reverse the mean shock velocity. Away from this 
special manifold our result for the sharpness of the shock suggest that finite 
systems with lattice size of the order 10 can be well described by the domain 
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wall theory for first order boundary-induced phase transitions, [HI El, with 
limitations analogous to those obtained from the exact results of Ref . [30] . 

It is intriguing that the maximal sharpness appears at some specific value 
of the driving force or, equivalently, ratio of boundary fugacities. It would 
be interesting to investigate whether such a field-induced sharpening of the 
interface is a special property of lattice models or can appear also in contin- 
uum systems such as the recently studied mass transfer models |3El dlj ■ It 
is also an open problem whether there can be an analogous reduction of the 
shock dynamics to a random walk problem in exclusion processes where the 
stationary distribution does not factorize [31^ I3T?]. 
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